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INVESTIGATION OF TC-1 FLIGHT FAILURE 


USING POWER SPECTRAL ANALYSIS 
by Carl F. Lorenzo 

Lewis Research Center 
Cleveland, Ohio 44135 

SUMMARY 

During the Titan Centaur 1 test flight a failure involving at least one of 
the Centaur propellant boost pumps occurred. Also, neither of the boost 
pump speed instruments indicated pump rotation. 

Accelerometer data from the Titan Centaur 1 flight failure were analyzed 
using power spectral density methods to determine boost pump speed during 
attempted starts of the Centaur. The technique was demonstrated on a ref- 
erence flight . 

The hydrogen boost pump speed transient was determined for the TC-1 
flight. Other trends are seen in the data. However, these are not believed 
to be the oxygen boost pump. Discussion of data enhancement techniques is 
also presented. 


INTRODUCTION 

During the Titan Centaur 1 test flight a failure of the Centaur vehicle 
occurred. Grossly, this failure involved at least one of the Centaur propellant 
boost pumps. Also, during this flight neither of the boost pump speed mea- 
surement instruments gave any indication of boost pump rotation. 

In view of the suspected boost pump failure, it was desired to use other 
telemetered data to determine as far as possible the boost pump speed tran- 
sients during the attempted starts of the Centaur stage. 

It was felt that the use of power spectral density methods would allow 
maximum use of existing data in terms of determination of boost pump oper- 
ation. Further, it was felt that the boost pump speed transient information 
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might be helpful in determining the suspected failure mode of the Centaur 
boost pump(s). 


BASIC SCHEME 


The basic philosophy was to apply the power spectral approach to an 
earlier successful reference flight in which the boost pump speed indicators 
worked properly. Then the results obtained using power spectral density 
methods would be compared with those indicated directly by the instruments, 
thereby validating the power spectral approach. Then, having this tool in 
hand, to apply it to the TC-1 flight data with greater confidence than would 
be had by directly applying the techniques. The reference flight selected 
was the Atlas/Centaur-31 flight. 

For purposes of this study, the following quantities are of interest: 

The autocorrelation function R (t) is defined by 

A 


R (t) = lim — 

A rp 

T — oc 



x(t)x(t + r)dt 


The power spectral density function G x (f) is defined by 


f" 


G x (f) = 2 I R^rJe^^dT = 4 j R x (t)cos 27rfr dr 




'0 


This can also be expressed as 


G*«> - 


lim lim — — 
Af — 0 T — <* ( Af ) T , 



(t,f, Af)dt 


where x(t, f, Af) is that portion of x(t) in the frequency range of f to f + Af 
(as obtained by passing xft) through a sharp cutoff filter of bandwidth Af at 
frequency f). 
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For more details and applications of these functions, see ref. 1 or 
similar texts. The power spectral density function G (f) will have peaks at 
multiples of one per revolution frequency. The basic data display then will 
be the so-called 3D or Waterfall plots in which individual plots of power 
spectral density versus frequency are stacked for successive time slices 
through the data to give a picture over a desired time period of the frequency 
content of the signal. 

In order to define the pump speed transients with sufficient resolution, 
the time slice required for each power spectral density plot was necessarily 
very short. This, therefore, dictates a correspondingly high frequency 
range for the PSD analysis. The analysis times used in this study were 1 
second slices which corresponded to a 500 Hz frequency range, and 1/2 sec- 
ond which corresponded to a 1000 Hz range (display truncated to 750 Hz). 

The frequency range was considerably above the expected one per revolution 
of the boost pumps in order to use the harmonic structure to identify the 
boost pumps against the background noise. 

REFERENCE ANALYSIS - ATLAS CENTAUR-31 FLIGHT 

A variety of accelerometer locations were edited to determine those 
most likely to expose boost pump speeds for the Atlas Centaur 31 (AC -31) 
flight data. Considered first were accelerometer locations common to those 
of the TC-1 flight at the forward end of the Centaur. 

That the distant (forward end) accelerometers were not useful is prob- 
ably due to the fact that the structural transmission path from the thrust 
section to the front end of the Centaur would only allow information to be 
passed from the pumps at a few frequencies associated with the resonances 
in the structural transfer functions between these two points (namely the 
transmissibility windows). Several pressures were measured in both the 
AC-31 flight and the TC-1 flight; however, the frequency bandwidth of the 
measurement was so low as to preclude their use for this study. 

The most likely results for identifying the speed transient on AC -31 
were obtained on accelerometers located in the thrust sections of the Centaur 
stage. Accelerometer locations are shown in fig. 1. The data presented 
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herein are taken from a single accelerometer CA880 in the thrust section. 
Unfortunately, there was no accelerometer at this location during the TC-1 
flight; however, the data quality of this accelerometer appeared superior to 
the remaining choices for this application. Therefore, commonality was 
sacrificed. 

The primary results of the power spectral analysis of accelerometer 
CA880 for the AC -31 flight are shown in figs. 2 and 3. The frequency range 
is 0 to 500 Hz in all cases and reference times are marked beside the photo- 
graphs. The photographs also provide some time overlap. 

The frequency concentration associated with the boost pumps are readily 
discerned in the photographs. Also, the harmonic relationship between the 
frequencies is easily seen.* The constant frequency peaks which are seen in 
the data are usually structural frequencies (perhaps local) which are constant 
over sustained periods. The b photographs of the series are marked to show 
the locus of the PSD peaks. The spectral peaks were read from these plots 
and tabulated. The lowest frequency mode seen in the figure has been inter- 
pretted as the second harmonic of the liquid oxygen boost pump. The next 
higher frequency (not harmonically related) has been interpretted as the 
second harmonic of the liquid hydrogen boost pump. The true speed lines as 
indicated by the speed measuring instrument on the pumps were then deter- 
mined. The results of these two approaches were then plotted as functions 
of time (fig. 4 and fig. 5). 

It is particularly interesting in this flight that the one per revolution 
frequency of the hydrogen boost pump at steady state occurred at the same 
frequency as the two per revolution frequency of the oxygen boost pump. 

This may indicate some difficulty at separating pumps for the TC-1 flight or 
in identifying a single pump. 

Several observations are made concerning the analysis of data for the 
AC -31 flight. First of all, in observing the power spectral density plots 
(figs. 2 and 3) directly, it is noticed that the low frequency end contains a 
significant amount of noise. This low frequency content is apparently engine 

*On viewing these photographs, often the patterns may be more easily 
recognized by viewing from the end of the pap#*- at small incidence angles. 



5 


noise propagated by the structure from the still attached Atlas booster. This 
noise tends to mask some of the lower frequency power concentrations gen- 
erated by the boost pumps rotation. 

The plots comparing measured speed and PSD implied boost pump speed 
(figs. 4 and 5) compare very well in shape, but differ somewhat in time. 

There are several effects which are probably influencing this time mismatch. 
One of these is accurately starting the power spectral analyzer for the 3D 
plots which is a manual operation based on a tape time signal This is 
probably the most important effect. The second effect is a question of sweep 
rate effect. That is, since the forcing functions (pump speeds) are varying 
in time and we are sensing an acceleration which is the result of some struc- 
tural transfer function, we have essentially the sweep rate effect which oc - 
curs in a vibration test; namely, the response appears to be delayed in time 
due to the time necessary for the structural transfer function to reach full 
value (sinusoidal steady state). 

The general results of figs. 4 and 5 appear to be a satisfactory matching 
when temporal effects are neglected. It was felt that the techniques could be 
satisfactorily applied to the TC-1 flight data. 

TITAN/CENTAUR 1 FLIGHT ANALYSIS 

For the Titan/Centaur 1 (TC-1) flight data, a variety of accelerometer 
locations were edited and studied. The accelerometer CA300 was con- 
sidered to be the best measurement in terms of boost pump activity standing 
above the background noise level. 

The spectral analysis of the (TC-1) flight is presented in figs. 6 through 
14, which covers the period prior to the first attempt to start the engine to 
the time after second start attempt. Again, the b part of the figures in- 
dicates the interpretation of boost pump power spectral concentrations and 
the associated harmonics. The data of fig. 11 was not used in the speed 
transient plots but is included here for completeness. 

The analysis frequency range of 0 - 1000 Hz was used. This corresponds 
to a 1/2 second analysis time - the time for each spectral scan in the plots. 
Only the frequency range of 0 - 750 Hz was displayed in the figures. The 
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lowest frequency interpretation lines have been interpretted as the second 
or fourth harmonics of the implied pump speed transients. 

The boost pump start signal occurred at 13:55:18. 3 for the first start 
attempt. Stage separation occurs at 13:55:56.7. 

It will be noted that the background noise level changes significantly 
after separation when the Centaur stage does not have the noise inputs con- 
tributed by the Titan engine system. 

Only a single pump is readily discernible in the plots that are shown. 
(That is, the curves are in harmonic relationship. ) The determination of 
which pump is being observed in the spectral data can be based on several 
discriminators. These are 1) Harmonic structure; strength of the various 
harmonics. 2) The rotor time constants as compared to previous flight; and 
3) Use of the hydraulic performance of the pumps in the system (i. e. , steady 
state pressure measurements). 

Some attempt was made to use the harmonic structures as a discrimin- 
ator, but the comparison was difficult and inconclusive. 

Based on previous flights of the Centaur, it was known that the time 
constant of the liquid oxygen boost pump was about 10 seconds, and that of 
the liquid hydrogen boost pump was 14 seconds. The time constant ab- 
stracted from the data for this flight was about 14. 5 seconds. In addition, 
the frequencies associated with the fundamental in this data matched very 
well with the expected one per revolution frequencies of the hydrogen boost 
pump. This result, together with the steady state hydraulic data allowed the 
conclusion that the pump observed in the spectral plots was the hydrogen 
boost pump. On that basis, the results of these transients interpreted as 
pump speed has been plotted in figs. 15 and 16 against time in flight for both 
the first and second start attempts for the Centaur engine system in the TC-1 
flight. Also indicated on these plots are various important times in the flight 
in terms of the engine start sequence . 

It will be noted that because of the masking noise of the Titan engines 
the speed transient cannot be determined at low speeds (frequencies). Also, 
the steady state speed is a little low compared to the AC -31 flight. That is 
7500 rpm compared to 8500 rpm. 
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LIQUID OXYGEN BOOST PUMP 

Several patterns were seen in the spectral plots which were considered 
to be possible indicators of liquid oxygen boost pump rotation activity. In 
general, any plots in which at least two harmonic lines were not indicated 
was rejected as possible lox boost pump indicator. 

Two plots which might represent lox boost pump rotation are shown in 
figs. 17 (fig. 7 repeated) and 18. In both figures the patterns indicated are 
at the limits of optical perception, and may not be valid. They are presented 
here only to indicate possible behavior of this lox boost pump. 

In fig. 17 the frequency relationship between the two patterns indicated 
is about 7 to 4. The maximum speed based on the 320 Hz level of the fourth 
harmonic is 4800 rpm. It is extremely unlikely that this is the lox pump, 
since the pressure data does not validate such a speed. In fig. 18 the har- 
monics are in a 2 to 3 relationship; here the maximum frequency is again 
320 Hz, which for the third harmonic implies 106 Hz fundamental which re- 
lates to a 6360 rpm rotational speed. 

It is not likely that either of these patterns is representative of the liq- 
uid oxygen boost pump, since the speeds are much higher than would be ex- 
pected from hydraulic performance. It is believed that these are probably 
the result of some rotating machinery in the still attached Titan vehicle. 
Again, the patterns are extremely weak and are at the limits of human per- 
ception. It is felt that some technique to enhance these results, to separate 
out the speed changes from background noise from the Titan engine is de- 
sired. How this might be done is explored in the Appendix. 

CONCLUDING REMARKS 

It has been demonstrated that the use of power spectral methods applied 
to accelerometer flight data can be used to indicate pump speed by testing 
the technique against known flight results for the AC-31 flight. 

The transient of hydrogen boost pump speed versus time for both at- 
tempted engine starts for the TC-1 flight has been established. The results 
appear to be quite good over most of the range, however, they suffer slightly 
at low amplitudes due to background noise. 
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There are patterns in the data of other frequency concentrations which 
may be other rotating machinery in the vehicle; however, it is not believed 
that these are representative of the oxygen boost pump. In addition, these 
patterns are very close to the noise level and are quite difficult to establish 
as valid. 

Some enhancement technique could be applied to these data to remove 
the effects of Titan engine noise; to either expose the oxygen boost pump 
speed transient (if it exists) and/or the hydrogen boost pump transient at 
the low speed condition. Such enhancement procedures have been outlined 
in the Appendix. 
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APPENDIX 

Some consideration was given to enhancement techniques which could be 
applied to the TC-1 data. These techniques were not implemented due to 
hardware problems in transferring the analyzed data from the spectrum an- 
alyzer to computers where the enhancement procedure could be performed. 
However, the techniques which were to have been applied are presented here 
for their potential value for future problems. 

The enhancement problem can be under sto*" d by reference to fig. 19. 

Here it can be seen that the accelerometer being analyzed is receiving in- 
formation from the two boost pumps .long with the undesired information 
from the Titan stage 11 which is still attached during the early part of the 
boost pump start transient. The information from the Titan consists of 
both engine and pump noises as transmitted (colored) by the structural 
transfer functions. 

However, since information is available prior to boost pump start (with 
only the Titan stage n information arriving at the accelerometer), it should 
be possible to form a reference spectrum during this period which could then 
be removed from the data during the boost pump start transient. This process 
should then yield an enhanced result. 

Several techniques to accomplish this are outlined in the paragraphs 
which follow. In the following material: 

is the instantaneous i th (time) unenhanced spectrum 
Ggj is the enhanced spectrum at the i th time 
z id 

Gj^ is a reference spectrum to be defined for each approach 
K is a constant 


TECHNIQUE 1 


For this technique the enhanced spectra is given by the equation: 
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Jg i - KG r for G t - KG r > 0 
Vo for G { - KG r < 0 


(Al) 


where K is greater than zero and will normally be chosen to be less than 
one. For this case, G R is a reference spectrum formed prior to boost 
pump start when only Titan engine noise is being generated. The reference 
spectrum G R could be an average of several spectra taken during that per- 
iod or a single representative spectrum at the choice of the analyst. 

The basis of this procedure is to eliminate the effect of the Titan engine 
noise from the boost pump start transient. Namely, prior to boost pump 
start an average spectrum called G R would be generated as 

n 

G R = IY G (A2) 

n 

i=l 


It should be noted that this procedure is somewhat artificial in that KG R 
can be greater than Gj which would yield Ggj < 0 without limiting. With 
limiting, this can produce flat spots (zero) in the enhanced spectr~ . 


TECHNIQUE 2 

For this procedure the enhanced spectra Gg^ are formed as 

G - °i 

Ei i + kg r 


(A3) 


where G^ and G R are formed as in technique 1. In this case it is clear 
that the spectral value of G Ei can never be negative and can never be great- 
er than G^. This appears to be a much more appropriate manner of elimin- 
ating the unwanted noise. 

It will be noted that when G R * °» G Ei G i and when G R is large (for 
example, because of structural resonance coloring of Titan engine noise), 

Ggj is reduced in spectral power. This approach does not require the arti- 
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ficial limiting used in the previous technique. 

As a guideline in selecting the value of K, define a reduction ratio R as 


1 + KG R MAX 
1 + KG R MIN 


(A 4) 


then K is determined to be 

K = (A 5) 

G R MAX “ RG R MIN 

The reduction ratio R represents the relative reduction of the peak in the 
reference spectrum relative to the minimum of that spectrum. 

It is also noted that the enhancement procedure Gg i = (Gj/G^) is a 
special case of Eq. (A3) where K is very large. 

TECHNIQUE 3 


In this case, 


G Ei " G i “ KG R 


(A6) 


however, G^ is now defined as 


Gp> (G: + 

Ri 2n i-n 


4 G i-2 4 G i-1 4 G i 4 1 4 


+ G 


i+n) 


(A7) 


The objective here is to form a reference spectrum which is representative 
of the nonchanging part of the 3D display. Removal of this reference spec- 
trum from G^ as indicated in Eq. (A6) should expose those parts of the boost 
pump speed transient which are changing most rapidly, namely the initial 
part of the transient. With this approach, when the boost pump speed reaches 
steady state it would be eliminated from the result. In essence, this tech- 
nique would expose the non- stationary part of the data. 

A variant to this approach would be to weigh the spectra composing G Ri 
for example 
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G Ri " ^ (A i-n G i-n + * * * * A i-2 G i-2 + A i-l G i-l + A i G i + * * * ’ A i4n G i40 ) 


or 


140 


G . = ... — (A. G. . . . . Ai jK'Gi ) = — 
ri 2 n - 1 1-n 1_n 1+11 * +n 2 


TECHNIQUE 4 


— V AG 

•*> z. " 


, (A8) 


J=l-n 


The above reference spectra (technique 3) can also be applied to tech- 
nique 2 where 



G i 


1 + KG 


R 


(A9) 


Here again only the chaining part of the spectra will be exposed, and 
the possibility of Gg < 0 does not exist. 
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